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ABSTRACT 


Application of active control to separated flow on the 
RC ( 6 ) -08 airfoil at high angle of attack by localized surface 
heating is numerically simulated by integrating the compressible 
two-dimensional nonlinear Navier-Stokes equations solver. Active 
control is simulated by local modification of the temperature 
boundary condition over a narrow strip on the upper surface of 
the airfoil. Both mean and perturbed profiles are favorably 
altered when excited with the same natural frequency of the shear 
layer by moderate surface heating for both laminar and turbulent 
separation. The shear layer is found to be very sensitive to 
localized surface heating in the vicinity of the separation 
point. The excitation field at the surface sufficiently altered 
both the local as well as the global circulation to cause a 
significant increase in lift and reduction in drag. 



INTRODUCTION 


The first known effort on the use of boundary layer control 
for reducing drag in the U.S.A. was made by NACA at the Langley 
Laboratory in the Low Turbulence Pressure Tunnel ( LTPT ) in 1939. 
In 1941, NACA conducted what is believed to be the first laminar 
suction experiment. This was reported in an unpublished paper by 
Maskell and described by Flatt 1 in 1961. This effort was 
apparently triggered by the results of contemporary German work, 
later reported by Holstein*. 

In recent years, the application of boundary layer control 
has been divided into high lift and low drag. The first deals 
with separation control and the second with laminar flow control 
and turbulent manipulation* - 1 8 . 

Design constraints between high speed cruise and low speed 
take-off and landing has forced the designers of fixed wing 
aircraft to incorporate some form of high lift device to improve 
lift at low speed. The improvements have been accommodated 
within the configuration of the trailing edge flap and leading 
edge slat. 

In rotary wing aircraft additional constraints exist 18 . In 
flight, the rotor blade is exposed cyclically to different flow 
about the azimuthal position in the rotational disk. The blade 
advancing into the oncoming airstream experiences high dynamic 
pressure and needs only small angles of attack to produce ade- 
quate lift to maintain roll trim. At the opposite side of the 
disk, the retreating blade experiences low dynamic pressure and 
therefore operates at high angles of attack to develop an ade- 
quate lift. When the airfoil reaches high angles of attack, the 
generated flow field can be characterized by massive unsteady 
separation originating near the leading edge and by large scale 
vortical structures. Within this complex flow, one encounters 
the phenomenon of dynamic stall, a behavior dominated by the 
dynamic stall vortex. 

Several investigators have developed analytical models for 
rotor flow 17 . Most cases provide a solution to the full 
potential equation. An important design consideration for the 
rotor and the blades is to operate at high lift conditions. 
Ideally, in this situation, the boundary layer must not separate. 
However, if the angle of attack is increased sufficiently, 
separation is unavoidable. Using flow control techniques, a 
number of benefits can be envisioned to control separation. One 
of these is to improve lift and stall characteristics. 
Applications of active control both in the laboratory 8 ' 7 and in 
numerical experiments 8-10 have been performed successfully. 



In this paper, the steady and unsteady flow about the 
RC( 6 ) -08 rotor airfoil relevant to the issue of dynamic stall and 
its origins and the method of controlling it by active means is 
investigated. In the simulation, the full two-dimensional 
compressible Navier-Stokes equations solver developed by Rumsey 
and Anderson and Rumsey and Thomas 1 • ' 1 * is used and the control 
is by surface heating. Their work successfully predicted lift 
and stall angle for airfoils. The flow being investigated is 
characterized by the shedding of a strong vortex from the leading 
edge. The viscous layer is of the order of the airfoil chord 
during the vortex shedding process. This class of flow problems 
may properly be simulated only by use of the full Navier-Stokes 
equations . 

The results presented herein contains both steady and 
unsteady data including lift and drag coefficients at Reynolds 
numbers (based on chord) Re of 10“ and 3x10* for laminar simu- 
lation and 5.2x10* for turbulent simulation. The objective is 
to investigate the active control of disturbances in the boundary 
layer by localized surface heating. The concern is with the 
nature of the interaction of the time dependent disturbance in a 
separated boundary layer and the ability to actively perturb the 
layer to reduce and prevent separation. It is shown that control 
of separated flow can be successfully achieved using localized 
active surface heating. 


NUMERICAL METHOD 

The solution of the compressible two-dimensional 
Navier-Stokes equations is used to compute the laminar flow and a 
suitable turbulent modeling is used to compute the turbulent flow 
about the RC(6)-08 airfoil. The method is an upwind-biased 
implicit approximate factorization Navier-Stokes algorithm 
described by Rumsey and Anderson 1 * and Rumsey and Thomas. 1 * The 
complete form of the two-dimensional compressible Navier-Stokes 
equations (at low Reynolds number for laminar flow and high 
Reynolds number for turbulent flow) are evaluated. The equations 
are written in generalized coordinates and conservation form: 
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where n is the coordinate along the body, C is the coordinate 
normal to the body, p is the density, u and v are the x and y com- 
ponents of velocities, and e is the total energy. 

The viscous terms on the right hand side are given by 
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All viscous terms are centrally differenced and implicit cross- 
derivative terms are neglected in the formulation. This 
algorithm is accurate to the first order in time and second in 
space. Implicit spatial derivatives of the convection and pres- 
sure terms are first order accurate. For turbulent calculations, 
the two layer algebraic eddy viscosity model of Baldwin and 
Lomax 20 is employed. 
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Boundary conditions are applied explicitly. No slip, 
adiabatic wall conditions, as well as zero pressure conditions 
are applied on the body: 


u = v = 
3p 3a 2 

3C = 3C 


0 

= 0 


( 6 ) 


where a 2 is proportional to the temperature and a is the 
speed of sound. In the farfield, the subsonic freestream 
boundary conditions are determined through a characteristic 
analysis normal to the body 1 '. 


RESULTS AND DISCUSSION 

The numerical results are obtained for the RC(6)-08 airfoil 
at a freestream Mach number M„ of 0.3 and a chord Reynolds number 
Re of 10* and 3x10* for the laminar Navier-Stokes code. in 
addition, the turbulent version of the code is evaluated at 
Mo„=0.35 and Re=5.2xl0‘. The first part demonstrates the laminar 
flow modeling and the active control procedures. The second part 
demonstrates the turbulent flow modeling, a comparison with wind 
tunnel data and the active control procedures. 

Active control is modeled by surface heating which modifies 
the temperature boundary condition over a small portion of the 
airfoil surface. The width of the control strip (h) is about 10 
percent of the wavelength of the flow disturbance at the leading 
edge. The control temperature is modified according to the 
formula 


T(t) 


ref 



+ A + B 2 sin f n t 


(7) 


In (7), T w is the temperature of the wall (520 R° ) and T ^ is 

the reference temperature (460 R°) where f, t, and n are the 
frequency, time, and an integer, respectively. The function 
models a D.C. and an A.C. current. The parameters of the heater 
for laminar flow are A=0.1 and B=1.22 for Re=l0 5 , and A=0.1 and 
B=0 . 87 for Re=3xlO*. For turbulent flow at Re=5.2xl0‘, A=0.1 
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and B=0.61 are chosen. The maximum temperature of the heating 
strip is approximately l . 5 times the unheated absolute wall 
temperature. This temperature is lower than the temperature 
reported in the experiment* using tungston wire. 

Accuracy of the computer solution and physical realism 
(i.e., flow patterns consistent with experimental observations) 
of the flow about the airfoil have been established by Rumsey and 
Thomas, including grid density and extent studies. In the pre- 
sent study, the grid refinement extending 10 chord lengths indi- 
cated that for laminar flow modeling a C-grid of 251x100 was 
sufficient while with turbulent flow modeling (up to 15 degrees 
angle of attack) a C-grid of 301x100 and 601x100 was used on a 
Cyber-205 computer and a Cray-2 computer respectively. The time 
step for the 301x100 grid was 0.001 and for the 601x100 grid 
was 0.0005. 


LAMINAR FLOW AND CONTROL 

The airfoil model is a helicopter blade section as depicted 
in Fig. 1. For laminar flow, consider first the freestream Mach 
number Mo, = 0.3 and an angle of attack a=10°. The separation 
region starts at the leading edge and continues downstream over 
the chord and into the wake region exerting large scale periodic 
oscillation. 

Active control of the separation is accomplished by 
modulating the natural shedding frequency with surface heating. 
Tables 1 and 2 and Figs . 2 to 7 show the results of the 
uncontrolled and controlled cases. The steady state, the root 
mean square (RMS), and instantaneous lift and drag 

coefficients Cl, Cq, */Cl 2 , -/Cp 2 , CL(t), and Cpft) show an 
increase in lift and reduction in drag for the control runs. 
Surface heating induces both local as well as global circulation 
about the upper surface. The global field is characterized by 
the vortices being shed periodically from the separation point at 
the leading edge. These vortices trigger large scale circulatory 
motion with clockwise rotation upstream of the trailing edge and 
with counterclockwise downstream. The separation point at the 
leading edge displays characteristic instability in the form of a 
sequence of discrete vortices consistent with experimental 
observations 7 . Using the conventional definition for Strouhal 
number S(f)=(f c sina)/U where U is the free stream velocity, 

the fundamental mode S(f]_) is 0.13 and the first harmonic mode 
S ( f 2 ) is 0.26 (Table 1). These results are in close agreement 
with that for a Joukowski airfoil 19 . 

The normalized perturbation velocity at half chord is shown 
in figures 4 and 5. The maximum perturbation occurs at 2 percent 
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of the chord in the vertical plane and its period is the same as 
the global period. The RMS increase due to surface heating is 
the difference between 0.307 and 0.449 and between 0.288 and 
0.379 for the two Reynolds numbers. The strong coupling between 
surface heating and separated flow causes the shear layer to 
resonate through its harmonics inducing additional vortical 
motion that remains close to the airfoil. This feature contrib- 
utes to an increase in lift. The normalized mean velocity pro- 
files (Fig. 3) show a definite reduction in separation affecting 
all the upper surface. Furthermore, the pressure distribution 
(Fig. 6) shows the region in which control has been most effec- 
tive over the airfoil. The functional form of the heat input is 
plotted in Fig. 7. The period matches that of vortex shedding. 
The summarized results in Table 2 demonstrate that localized sur- 
face heating has obvious attraction as a means for actively 
exciting the region of flow separation. Results confirm that an 
externally imposed temperature field increases the perturbation 
(as shown by the higher RMS value) within the separated layer and 
thus increases the compactness of the flow and reduces the 
separation. 


TURBULENT FLOW AND CONTROL 

This section is divided into two parts. The first part 
discusses the comparison of analysis with experiment for pressure 
distributions, lift coefficients and drag coefficients at 
different angles of attack at a Reynolds number of 5.2x10* and a 
Mach number of 0.35. In the second part, the effect of active 
control by surface heating is presented and discussed. 


1. Comparison with Experiment 

The comparison of lift and drag coefficients and pressure 
distributions from the Navier-Stokes code with wind tunnel 
measurements (unpublished data) is shown in Figs. 8 and 9. The 
experiment provides a broad range of conditions up to an angle of 
attack a=12.46 0 . The major difference between the experiment and 
the numerical result is in the pressure distribution near the 
leading edge at high angles of attack. The lift coefficients 
increase linearly up to a~8° and beyond that the coefficients 
exhibit a behavior resembling that of a stall of mixed type. The 
comparison of Cl and Cq are accurately predicted about the mid- 
dle range of the angles of attack. The difficulty for the larger 
a condition where flow separates is in the resolution of the sur- 
face pressure since the perturbation field is a significant por- 
tion of the total signal. In addition, the experimental values 
of the maximum lift coefficients of an airfoil measured in the 
Langley 6- by 28-Inch Transonic Tunnel are known to be low (on 
the order of 0.10) due to the influence of the tunnel sidewall 
boundary layers. 21 The resolution at higher angles of attack in 
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the presence of large separation needs further study both from 
the experimental and numerical points of view. From a numerical 
point of view, further grid resolution is needed at higher a. 

For the experiment, a real-time resolution from D.C. to higher 
frequencies is also needed. 


2. Active Control 

Results from the controlled and uncontrolled runs are shown 
in Figs. 10 through 16. The width of the control strip (h) is 
2 percent of the airfoil chord and the strip is located near 
the leading edge Fig. 1. The control temperature T(t) reaches a 
peak at about 1.5 times the local steady ambient value. At the 
highest angle of attack, a=15°, the computation is made with two 
different grid sizes 301x100, and 601x100, since it was found 
that grid refinement is needed for a>13° because of the 
separation which dominates over the upper surface. 

The perturbation velocity at the half chord point at a=l5° 
is lower closer to the surface and it increases with distance 
from the surface, an indication that the flow is separated over 
part of the layer. Surface heating activated near the leading 
edge augments the circulation over the upper surface thus 
reducing the separation significantly. This is observed by the 
increase in velocity perturbation. The improved circulation is 
also apparent by the mean velocity profile at three stations on 
the upper surface where the controlled profile has significantly 
reduced the region of separation as compared with the 
uncontrolled profile Fig. 11. The velocity perturbation is 
characterized by discrete harmonically related peaks with 
nondimensional frequency fc/a^O . 145 ; 0.290; 0.445; 0.630; 0.760 

and 0.950 at y/c=10~* and above while close to the surface only 
the odd harmonics are excited. The spectra of the velocity for 
the control case exceeds the amplitude of the uncontrolled one at 
all frequencies, especially with increased distance from the 
surface. For the control run, frequencies tuned to the naturally 
shedding frequencies are used to enhance receptivity. When the 
power spectral density of the longitudinal velocity for the 
controlled case is compared with the corresponding uncontrolled 
case, the effect of surface heating is clearly noticed. The 
coupling between surface heating and separated flow causes the 
shear layer to resonate through its harmonics, thus inducing 
additional vortical motion as observed in the experiments. 6 ' 7 

The effect of control in the variation of C£ with a and 
with C£ is to provide a significant increase in lift and 
reduction in drag. Figs. 14 and 15. The effect of grid size is 
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also noticed for the higher a indicating an increase in grid 
resolution is needed with the increases in angle of attack. 

The pressure coefficient time history is given in Fig. 16. 

The amplitude at three different nondimens ional times 
( ta^/c ) =10 . 00 , 10.75, and 11.50, responds globally to the 

variation of the oscillating pressure about the airfoil. Regions 
of high pressure for the control run are maintained over the upper 
surface indicating a smaller separation than the uncontrolled run 
with lower pressure. This can also be noticed from the lift and 
drag coefficient time history Fig. 2 at lower Reynolds number 
illustrating that the uncontrolled runs maintain lower 
fluctuating lift and higher fluctuating drag over the cycle. 

Control of the unsteady turbulent separation is not yet 
well understood in a fundamental sense. The numerical analysis 
of the flow field continues to search for an efficient code that 
can be used over broader ranges of practical situations. 

The steady state compares with the experiment 
satisfactorily. The RC(6)-08 airfoil has good performance 
due to the extended unseparated region over the upper surface up 
to a~8 ° . It produces substantial high lift and low drag (at 
about 5 million Reynolds number, the upper limit of the present 
investigation). Active separation control works well on this 
airfoil since it delays massive separation and stall. 


CONCLUDING REMARKS 

The objective to show performance enhancement of an airfoil 
at high angle of attack was demonstrated at low and high Reynolds 
numbers by means of actively controlling the separation using 
localized surface heating. The results illustrate that the 
effect of control of separated flow is very much related to the 
wave form, the amplitude, and the bandwidth of the temperature 
function for a given pressure gradient.* The temperature 
function of the surface heating can restructure the lift and the 
drag of the airfoil profoundly. When both amplitude and temporal 
behavior are chosen for high receptivity, improved performance 
beyond the design limit can result. Previously, active surface 
heating has been shown experimentally to trigger harmonics and 
transition. Furthermore, it was observed that forced excitations 
induce instability and resonance and increase the entrainment in 
the early part of the shear layer, and as a result, the separated 
region was reduced. Further experimental and numerical work is 
needed to study three-dimensional effects on rectangular and 
nonrectangular configurations. 

There are many other practical applications in which some 
form of coupling — thermal, acoustic, or vortical exist to 
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enhance circulation and performance. This study has proven heat 
addition to be a very effective method for actively perturbing 
the separated shear layer resulting in a significant increase in 
lift and reduction in drag. 
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Figure 16 . Inscenceneous total pressure coefficient. 
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